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^ '■ Abstract 

I I Substantial collective flow is observed in collisions between large nuclei at RHIC (Relativistic 

^ I Heavy Ion Collider) as evidenced by single-particle transverse momentum distributions and by 



c 



azimuthal correlations among the produced particles. The data are well-reproduced by perfect 



^ , fluid dynamics. A calculation of the dimensionless ratio of shear viscosity rj to entropy density s by 

C^ I Kovtun, Son and Starinets within anti-de Sitter space/conformal field theory yields ry/s = h/AnkB 

which has been conjectured to be a lower bound for any physical system. Motivated by these 
results, we show that the transition from hadrons to quarks and gluons has behavior similar to 
helium, nitrogen, and water at and near their phase transitions in the ratio rj/s. We suggest that 
experimental measurements can pinpoint the location of this transition or rapid crossover in QCD. 

PACS numbers: 12.38.Mh, 25.75.-q, 25.75.Nq, 51.20.+d 



One of the amazing experimental discoveries of measurements on gold-gold collisions at 
RHIC (Relativistic Heavy Ion Collider) at Brookhaven National Laboratory is the surprising 
amount of collective flow exhibited by the outgoing hadrons. Collective flow is evidenced in 
both the single-particle transverse momentum distribution [Ij , commonly referred to as radial 
flow, and in the asymmetric azimuthal distribution around the beam axis J3|, quantified by 
the functions vi{y,pT),V2{y,PT), ••• in the expansion 
d^N 1 (PN 

where y is the rapidity and p^ is the transverse momentum. The function v^iy = 0,Pt), in 
particular, was expected to be much smaller at RHIC than it is at the lower energies of the 
SPS (Super Proton Synchrotron) at CERN |3], but in fact it is about twice as large. Various 
theoretical calculations jjj support the notion that collective flow is mostly generated early 
in the nucleus-nucleus collision, and is present at the partonic level before partons coalesce 
or fragment into hadrons. Theoretical calculations including only two-body interactions 
between partons cannot generate sufficient flow to explain the observations unless partonic 
cross sections are artificially enhanced by more than an order of magnitude over perturbative 
QCD predictions J5|. This has emphasized that the quark-gluon matter created in these 
collisions is strongly interacting, unlike the type of weakly interacting quark-gluon plasma 
expected to occur at very high temperatures on the basis of asymptotic freedom 6]. On the 
other hand, lattice QCD calculations yield an equation of state that differs from an ideal gas 
only by about 10% once the temperature exceeds 1.5Tc, where T^ ~ 175 MeV is the critical 
or crossover temperature from quarks and gluons to hadrons [Tf • Furthermore, perfect fluid 
dynamics (with zero shear and bulk viscosities) reproduces the measurements of radial flow 
and V2 very well up to transverse momenta of order 1.5 GeV/c jq|. 

An amazing theoretical discovery was made by Kovtun, Son and Starinets p| , who showed 
that certain special fleld theories, special in the sense that they are dual to black branes in 
higher space-time dimensions, have the ratio rj/s = 1/Att (we use units with h = k-Q = c = 1) 
where rj is the shear viscosity and s is the entropy density. They conjectured that all 
substances have this value as a lower limit, and gave as examples helium, nitrogen, and 
water at pressures of 0.1 MPa, 10 MPa, and 100 MPa, respectively. Interesting enough, this 
bound is also obeyed by A/" = 4 supersymmetric SU{Nc) Yang-Mills theory in the large Nc 
limit |lo|. 



We are motivated by these discoveries to study what happens in QCD at finite temper- 
ature. The relatively good agreement between perfect fiuid calculations and experimental 
data for hadrons of low to medium transverse momentum at RHIC suggests that the vis- 
cosity is small; however, it cannot be zero. Indeed, the calculations within anti-de Sitter 
space/conformal field theory suggests that t] > s/{4tt). Our conclusion will be that suffi- 
ciently precise calculations and measurements should allow for a determination of the ratio 
rj/s as a function of temperature, and that this ratio can pinpoint the location of the phase 
transition or rapid crossover from hadronic to quark and gluon matter. This is a different 
method than trying to infer the equation of state of QCD in the form of pressure P as a 
function of temperature T or energy density e. 

The energy-momentum tensor density for a perfect fluid (which does not imply that the 
matter is non-interacting) is T^'^ = —Pg^'^ + wu^u^ . Here w = P + e = Ts is the local 
enthalpy density and u^ is the local flow velocity. Corrections to this expression are pro- 
portional to flrst derivatives of the local quantities whose coefficients are the shear viscosity 
r] and bulk viscosity C,. (Thermal conductivity is neither relevant nor defined when all net 
conserved charges, such as electric charge and baryon number, are zero.) Explicit expres- 
sions may be found in textbooks [ll|. Perfect fiuid dynamics applies when the viscosities are 
small, or when the gradients are small, or both. The dispersion relations for the transverse 
and longitudinal (pressure) parts of the momentum density are 

u + iDtk'^ = 
u^ -v^P + iDiujP = (2) 

where Dt = rf/w and Di = {^rj + C)/w are diffusion constants with the dimension of length 
and V is the speed of sound. Since w = Ts, and since usually the bulk viscosity is small 
compared to the shear viscosity, the dimensionless ratio of (shear) viscosity to entropy 
(disorder) rj/s is a good way to characterize the intrinsic ability of a substance to relax 
towards equilibrium independent of the actual physical conditions (gradients of pressure, 
energy density, etc.). It is also a good way to compare very different substances. 

In figures 1 through 3 we plot the ratio r]/s versus temperature at three fixed pressures, 
one of them being the critical pressure (meaning that the curve passes through the critical 
point) and the other ones being larger and smaller, for helium, nitrogen and water. The 
ratio was constructed with data obtained from the National Institute of Standards and 



Technology (NIST) [12] . (Care must be taken to absolutely normalize the entropy to zero at 
zero temperature: we did that using data from CODATA [13[.) The important observation 



14l Il5l | is that rj/ s has a minimum at the critical point where there is a cusp. At pressures 
below the critical pressure there is a discontinuity in r]/ s, and at pressures above it there 
is a broad smooth minimum. The simplest way to understand the general behavior was 
presented by Enskog, as explained in [16[. Shear viscosity represents the ability to transport 
momentum. In classical transport theory of gases rj/ s ~ Tli^eeV-, where /free is the mean free 
path and v is the mean speed. For a dilute gas the mean free path is large, /^ee ~ 1/na, with 
n the particle number density and a the cross section. Hence it is easy for a particle to carry 
momentum over great distances, leading to a large viscosity. (This is the usual paradox, that 
a nearly ideal classical gas has a divergent viscosity.) In a liquid there are strong correlations 
between neighboring atoms or molecules. A liquid is homogeneous on a mesoscopic scale, 
but on a microscopic scale it is a mixture of clusters and voids. The action of pushing on 
one atom is translated to the next one and so on until a whole row of atoms moves to fill a 
void, thereby transporting momentum over a relatively large distance and producing a large 
viscosity. Reducing the temperature at fixed pressure reduces the density of voids, thereby 
increasing the viscosity. The viscosity, normalized to the entropy, is observed to be the 
smallest at or near the critical temperature, corresponding to the most difficult condition 
to transport momentum. This is an empirical observation. The mean free path must lie 
somewhere between the dilute gas limit, 1/no", and the close-packing limit, 1/n^'^. For a 
massless gas of A^ bosonic degrees of freedom, with entropy density s = N{4:Tt'^/90)T^, the 
close-packed limit gives r]/s ~ 2/N^'^. 

How does this relate to hadrons and quark-gluon plasma? In the low energy chiral limit 
for pions the cross section is proportional to s/ f^, where s is the usual Mandelstam variable 
for invariant mass-squared and fj^ is the pion decay constant. The thermally averaged cross 
section is (a) oc T'^/ f^, which leads to rj/s oc {f.„/TY. Explicit calculation gives [12| 

^ = i^^ (3) 

s 167rT4 ^ ' 

Thus the ratio rj/s diverges as T ^ 0. At the other extreme lies quark-gluon plasma. The 
parton cross section behaves as a oc g'^/ s. A first estimate yields rj/s oc 1/ g'^. Asymptotic 
freedom at one loop order gives g^ oc l/ln(T/AT) where A^ is proportional to the scale 
parameter Aqcd of QCD. Therefore r]/ s is an increasing function of T in the quark-gluon 



phase. As a consequence, rj/ s must have a minimum. Based on atomic and molecular data, 
this minimum should lie at the critical temperature if there is one, otherwise at or near the 
rapid crossover temperature. 

The most accurate and detailed calculation of the viscosity in the low temperature hadron 
phase was performed in [l7|. The two-body interactions used went beyond the chiral ap- 
proximation, and included intermediate resonances such as the p-meson. The results are 
displayed in figure 4, both two flavors (no kaons) and three flavors (with kaons). The qual- 
itative behavior is the same as in eq. ©• The most accurate and detailed calculation of 

n 

the viscosity in the high temperature quark-gluon phase was performed in [1^. They used 
perturbative QCD to calculate the full leading-order expression, including summation of the 
Coulomb logarithms. For three flavors of massless quarks the result is 

1^ "■'' (4) 

a 9*lii(2.42/9) ^ ' 

We used this together with the two-loop renormalization group expression for the running 
coupling 

with At = 30 MeV, which approximately corresponds to using an energy scale of 2ttT and 
A.jjg = 200 MeV. The result is also plotted in figure 4. These results imply a minimum in 
the neighborhood of the expected value of Tc ^ 190 MeV. Whether there is a discontinuity 
or a smooth crossover cannot be decided since both calculations are unreliable near Tr,. 

It is interesting to ask what happens in the large Nc limit with g'^N^ held fixed [l^. In 
this limit, meson masses do not change very much but baryon masses scale proportional to 
Nc] therefore, baryons may be neglected in comparison to mesons due to the Boltzmann 
factor. Since the meson spectrum is essentially unchanged with increasing Nc, so is the 
Hagedorn temperature. The critical temperature to go from hadrons to quarks and gluons 
is very close to the Hagedorn temperature, so that Tc is not expected to change very much 
either. In the large N^ limit the meson-meson cross section scales as l/N^- According to 
our earlier discussion on the classical theory of gases, this implies that the ratio rj/s in the 
hadronic phase scales as N^. This general result is obeyed by (jH)) since it is known that 
/^ scales as Nc. The large Nc limit of the viscosity in the quark and gluon phase may be 



Q 



inferred from the calculations of [18| to be 

r/\ _ /l + 3.974r\ 69.2 

s/qgp~ V l + 1.75r j (^2;vj^ In (26/(^2 jv^(l + 0.5r))) ^^ 

where r = Nf/Nc- Thus the ratio rj/s has a finite large Nc limit in the quark and gluon 
phase. Therefore, we conclude that rj/s has a discontinuity proportional to A^^ if Nc — > cxd. 
This jump is in the opposite direction to that in figure 4. 

So far the only quantitative results for viscosity in lattice gauge theory have been reported 
by Nakamura and Sakai 20| for pure SU(3) without quarks. This bold effort obtained 
T]/s ~ 1/2 in the temperature range 1.6 < T/Tc < 2.2, albeit with uncertainties of order 



100%. Gelman, Shuryak and Zahed [2l| have modeled the dynamics of long wavelength 
modes of QCD at temperatures from Tc to 1.5Tc as a classical, nonrelativistic gas of massive 
quasi-particles with color charges. They obtained a ratio of rj/s ~ 0.34 in this temperature 
range. 

It ought be possible to extract numerical values of the viscosity in heavy ion collisions via 
scaling violations to perfect fluid flow predictions. One should perform a systematic beam 
energy and projectile/target mass scan from SPS energies to the top RHIC energy, and then 
on to the LHC Flow data, in the form of the functions fi,f25 •••; should be obtained and 
compared with the results of calculations based on relativistic viscous fluid dynamics. This 
program is analogous to what was accomplished at lower energies of 30 to 1000 MeV per 
nucleon beam energies in the lab frame. At those energies, scaling violations to perfect fluid 
dynamics were indeed observed J22]. It was possible to infer the compressibility of nuclear 
matter and the momentum-dependence of the nuclear optical potential via the transverse 
momentum distribution relative to the reaction plane 23| and via the balance between 
attractive and repulsive scattering J2J]. There is much to be learned about QCD at high 
energy densities. 
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(L. McLerran). 
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FIG. 1: The ratio r]/ s as a function of T for helium with s normahzed such that s{T = 0) = 0. The 
curves correspond to fixed pressures, one of them being the critical pressure, and the others being 
greater (1 MPa) and the other smaller (0.1 MPa). Below the critical pressure there is a jump in 
the ratio, and above the critical pressure there is only a broad minimum. They were constructed 
using data from NIST and CODATA. 
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FIG. 2: The ratio rj/s as a function of T for nitrogen with s normahzed such that s{T = 0) = 0. 
The curves correspond to fixed pressures, one of them being the critical pressure, and the others 
being greater (10 MPa) and the other smaher (0.1 MPa). Below the critical pressure there is a 
jump in the ratio, and above the critical pressure there is only a broad minimum. They were 
constructed using data from NIST and CODATA. The curves are plotted on logarithmic scale to 
make the behavior around the critical point more visible. 
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FIG. 3: The ratio r//s as a function of T for water with s normahzed such that s(T = 0) = 0. The 
curves correspond to fixed pressures, one of them being the critical pressure, and the others being 
greater (100 MPa) and the other smaller (10 MPa). Below the critical pressure there is a jump in 
the ratio, and above the critical pressure there is only a broad minimum. They were constructed 
using data from NIST and CODATA. 
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FIG. 4: The ratio r//s for tlie low temperature liadronic phase and for the high temperature quark- 
gluon phase. Neither calculation is very reliable in the vicinity of the critical or rapid crossover 
temperature. 
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